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Abstract 

A new method of employing an isospin chemical potential for QCD-like theories with different 
number of colors, number of fermion flavors, and in different fermion representations is proposed. 
The isospin chemical potential, which can be simulated on the lattice due to its positive definite 
determinant gives a means to probe both confining theories and IR conformal theories without 
adjusting the lattice spacing and size. As the quark mass is reduced, the isospin chemical potential 
provides an avenue to extract the chiral condensate in confining theories through the resulting pseu- 
doscalar condensate. For IR conformal theories, the mass anomalous dimension can be extracted 
in the conformal window through "finite density" scaling since the isospin chemical potential is 
coupled to a conserved current. In both of these approaches, the isospin chemical potential can 
be continuously varied for each ensemble at comparable costs while maintaining the hierarchy be- 
tween the lattice size and lattice spacing. In addition to exploring these methods, finite volume 
and lattice spacing effects are investigated. 

PACS numbers: 12.38.Gc, 12.39.Fe, 12.60.Nz 
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I. INTRODUCTION 



The study of strongly coupled beyond the Standard Model (BSM) theories through lattice 
gauge theory has seen a surge of activity within the last few years as a result of advancements 
in computational technology and the advent of the LHC In order to truly understand 
the non-perturbative dynamics that would lead to electroweak symmetry breaking in these 
technicolor theories [H E] (for reviews, see Ref. [31 0]), numerical lattice calculations are 
necessary. While not completely ruled out, scaled up theories of QCD with two flavors and 
three colors have been disfavored due to a larger than observed S-parameter [3 |6] , and a 
chiral condensate too small to account for the flavor hierarchy when extended technicolor 
and flavor changing neutral current bounds are considered [71 [8] . However, the dynamics 
of strongly coupled gauge theories with different number of flavors, colors, and fermion 
representations are still largely unknown. One proposed resolution for these experimental 
discrepancies is walking technicolor [5HT5]. which is thought to both yield an enhanced 
condensate and parity doubling. For fixed representation and colors, these scenarios are 
believed to set in for number of flavors just below the onset of the conformal window, 
after which spontaneous symmetry breaking no longer occurs. The majority of the lattice 
calculations to date focus on addressing these situations non-perturbatively. 

There have been a multitude of lattice studies involving asymptotically free QCD-like 
theories. The majority of calculations have explored three-color theories with fermions in 
the fundamental representation at various number of flavors [TBHSZ]- However, there have 
also been multiple studies of two color fundamental [2S], two color adjoint [2SHSS], and three 
color sextet [5BVHD] . While great progress has been made in each of these theories, much 
debate still remains largely due to the systematic uncertainties that result from discretization 
choices, non-zero quark masses, and finite volume effects. Also, unlike QCD, there is little 
phenomenology to guide these calculations. Thus, these calculations are heavily reliant on 
effective field theory (EFT) extrapolations or techniques for critical phenomena to attain 
information in the chiral limit. For these techniques to be valid, however, small fermion 
masses are required. 

In this work we explore the effect of introducing a new dimensionful parameter to the 
action, namely an isospin chemical potential, and how this parameter can help extract the 
chiral condensate in the chiral limit for confining theories and the mass anomalous dimension 
in the conformal window. The isospin chemical potential, which couples to the conserved 
isospin current, has been studied extensively in QCD [lllll2] and has been simulated on the 
lattice [131 - H6] . Isospin chemical potentials have also been explored as a probe for low-energy 
scattering parameters in meson-baryon systems [17| and two-point functions for mesons [l8] . 
In addition to not having a sign problem, it has several other beneficial properties. First, 
for a confining theory with chiral symmetry breaking, an isospin chemical potential leads to 
a pseudoscalar condensate, provided the chemical potential is larger than its critical value. 
The scale of this condensate is set by the confinement scale and can be related to the chiral 
condensate in certain limits. Second, unlike the fermion mass term, the isospin chemical 
potential is coupled to a conserved current, which results in the chemical potential not 
receiving any renormalization. This fact, along with a modest value for the mass, allows 
for one to extract the mass anomalous dimension via "finite density scaling" similar to the 
extraction of critical exponents from finite size scaling. A similar approach of exploiting 
the non-renormalization of the chemical potential to extract a critical index was performed 
in Ref. [IH]. Additionally, studies have been performed on both chemical potentials and 
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conformal phase transitions, along with free energy [50l |5T] . 

The organization of this paper is as follows. In Sec. [Tl| the isospin chemical potential is 
reviewed for both continuum QCD and lattice QCD. In Sec |III[ these results are generalized 
to confining theories and two methods are presented for extracting the chiral condensate 



in the chiral limit, namely, the pion mass and the pseudoscalar condensate. In Sec. IV 
universal finite density scaling functions, similar to the finite size scaling functions examined 
in Ref. [52H51] are derived for theories with an IR fixed point and methods for obtaining the 
mass anomalous dimension via finite density scaling are presented. In Sec. |V} an estimate of 
the lattice artifacts due small quark masses in the presence of an isospin chemical potential 
is given. 



II. ISOSPIN CHEMICAL POTENTIAL 



For each global symmetry/ conserved quantity that a system possesses, one can add a term 
to the Lagrangian which consists of a chemical potential coupled to the resulting conserved 
current. This term, which corresponds to the energy required to add an additional particle 
that carries one unit of the associated conserved charge, ultimately allows for exploration of 
the system at finite density. Physically, the system of the greatest interest is finite baryon 
density where a baryon chemical potential is coupled to the current of the conserved baryon 
number. However, due to to a complex fermion determinant in Euclidean space, the usual 
Monte Carlo lattice techniques cannot be applied to complex representation fermions^. In 
this work, we will explore a chemical potential coupled to a different conserved current, 
namely an isospin chemical potential. As long as the theory possesses an SU (2) subgroup 
of the flavor symmetry, this approach is applicable. 



A. Continuum QCD 



The Lagrangian for a two flavor QCD with an isospin chemical potential coupled to the 
third isospin matrix is given by 



C = ^ 



r3 



t{lp+ifino-) -M tP, (1) 



where ^ is the SU{Nc) covariant derivative, = diag(mq,mq), /i/ is the isospin chemical 
potential, and the fermion fields ip represent a pair of quarks in an isospin doublet. When 
/i/ = rrig = 0, the action has an SU{2)l ® SU{2)r chiral symmetry which breaks sponta- 
neously to the vector subgroup, SU{2)v- Turning on the chemical potential term, which 
behaves like the time component of a uniform gauge field, leads to the action being invariant 
under U{1)l ® U{1)r which spontaneously breaks to U{l)v (turning on the mass term gives 
an explicit breaking to this symmetry). This f/(l)y symmetry can be best understood as a 
circle in the plane perpendicular to the t^, term in the Lagrangian. 



Baryon chemical potentials are applicable to a host of theories that are in the real or pseudoreal represen- 
tation, such as two color theories, adjoint representation calculations, or SO{Nc) gauge theories 
due to a positive definite fermion determinant. While this could be a very interesting approach to many 
BSM calculations, we do not cover this scenario here. 
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One of the interesting outcomes to having an isospin chemical potential is the possibility 
of forming a pseudoscalar condensate [HI |32]. This condensate comes about when the 
the chemical potential is larger than a certain critical value and the U{l)v symmetry is 
spontaneously broken in the direction of the condensate. To understand this process, it is 
best to work with a small quark mass and small chemical potential in the chiral regime 
(namely, rriq ~ yU/ ^ A^, where is the chiral symmetry breaking scale). The leading 
order (LO) Lagrangian in chiral perturbation theory (xPT ) is given by [STf [58] 



tr(D^SD^St) + 2B tr(MtS + S+M) 



(2) 



where the normalization F ~ 130 MeV and M = diag(mg,mq). S is the matrix of the 
Nambu-Goldstone fields which can be parameterized as S = exp{ia n-r) = cos a+m-r sin a 
and B is related to the chiral condensate as the quark mass goes to zero. 



B = lim 



Ml 

F2 



(3) 



The isospin chemical potential enters in the covariant "gauge" derivative, D,^!! = d^H + 

3 

i[fijY^uo,^]- Expanding Eq. (2), one arrives at the potential 



tr( [r^ S] [r^ S"^]) - 25mgtr(S + S"^ 



(4) 



The potential as a function of a for several isospin chemical potential values is shown 
in Fig. [T] When hj < 2Bmq, the potential increases from its value at a = and vacuum 
remains unchanged. When fij > 2Bmq, the potential forms a new minimum and a condensed 
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FIG. 1: Chiral potential vs. a for values of /i/ below (blue), equal (red), and above (yellow) 2Bmq. 
When above the critical value, the potential develops a minimum at a non-zero value of a, which 
defines the vacuum expectation value in the condensed phase. 
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phase is favored. In terms of the parameter a, the minimum of the potential is given by^ 



2Bmg 

cos a = 2"-^. (oj 

As a result, the vacuum expectation value of the Nambu-Goldstone fields is given by 




> 2Bmq 



where n is the direction of the spontaneously symmetry. This new alignment will alter the 
chiral effective field theory and will consequently affect observables. This is accounted for in 
the xPT Lagrangian by expanding S about this new minimum, S = ^qS^o, where = ^o- 
This feature can be used to pick off low-energy coefficients, as proposed in Ref. ^Tj, and is 



explored for confining theories in Sec. Ill 



B. Lattice QCD 

When extending the continuum formulation of a field theory to a discretized simulation, 
there are often several subtleties that need to be understood. The first aspect is how the finite 
isospin density system maintains a positive definite determinant while the baryon chemical 
potential does not. One condition that can ensure this property for a theory with an even 
number of flavors is that the determinant of the action is real (multiplying an even number 
of real determinants is greater than or equal to zero). One relation that ensures a lattice 
action maintains a real determinant is if the discretized action, M, obeys the 75-hermiticity 
condition 

= -f^Mj^ — > det M = det MK (7) 

Consequently, the determinant of the product det M'^M = det M"^ det M = \ det is pos- 
itive definite. The inclusion of any chemical potential term of the form ^70^ will invalidate 
this relation. However, this is not the only condition to lead to a real determinant. One can 
include a flavor matrix and arrive at the relation 

= Ta^MTa^ — > detM = detM^ i e 1,2. (8) 

This relation is satisfied by isospin chemical potential term, ipr^'y^ip, ensuring that the 
determinant of the action is real. Note, that the baryon chemical potential coupled to ip'y^ijj 
will not satisfy this relation. 

Second, in a finite volume, spontaneous symmetry breaking without an explicit symmetry 
breaking source cannot occur, unless there are infinitely many degrees of freedom. The way 
to explore spontaneous symmetry breaking phases on the lattice is to include a small source 
of explicit symmetry breaking and extrapolate to when the source is zero. This is the case 
when exploring the chiral condensate since simulations often have an explicit chiral breaking 



^ In order to reduce confusion throughout this work, we leave the leading order vacuum pion mass as 2Braq 
(Gell-Mann-Oakes-Renner relation). The quantity will be reserved for the pion mass measured in the 
presence of the pion condensate. 
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mass term along with any chiral breaking discretization effects. Wliile tlie source for chiral 
symmetry breaking is inlierent to tlie system, tlie same cannot be said for a pseudoscalar 
condensate. Thus, in order to measure the condensate, a small breaking term is required 
to "tip" the condensate in a specific direction the the r^-r^ plane. The pseudoscalar source 
term added to the Lagrangian is of the form [17| 



A£ = zeV-.T y ^x, (9) 

which resembles a twisted mass term with a flavor coupling. In addition to dictating the 
direction of the pseudoscalar condensate, this term also is necessary to ensure the Dirac inver- 
sion does not become singular H5] . If one accounts for this term in the xPT Lagrangian, 
the minimum of the potential in Eq. Q W2\ 

25m„ Be . , 

cosa = ^ 5- cot a. (10) 

As long as e is small (taken to zero before rriq or /i), the continuum, infinite volume result 
is recovered. However, it is worth noting that if e is not sufficiently small, the cot a term 
diverges as the critical limit is approached from above, rendering the expansion invalid. The 
resulting vacuum expectation value (talking e — )■ first) leads to an altered chiral condensate 
and a pseudoscalar condensate given by 

(tjjtjj) = F^Bcosa, 

^(V'yTsV') = F'^Bsina. (11) 

Third, while discretization effects are less pronounced for the isospin chemical potential 
case than chiral symmetry breaking, care needs to be taken when performing tunings or 
interpreting residual mass terms (especially when taking rriq to zero). The system now has a 
pseudoscalar condensate and a chemical potential which lead to several of the alterations of 
the usual tuning methods (such as ensuring the pion mass vanish in the chiral limit) as will 



be described in Sec. |III[ Additionally, the isospin chemical potential term should correspond 
to the time component of the conserved vector current. Thus, in certain formulations, such 
as Domain Wall fermions, the used of the local vector current may not yield the desired 
results. However, this concern appears to be insignificant in lattice studies of small baryon 
chemical potentials [59j . 



III. CONFINING THEORIES 

Confining theories in the context of this paper are theories where dynamical symmetry 
breaking occurs and a confinement scale emerges. One class of theories similar to QCD, 
has scales that follow the hierarchy (including the isospin chemical potential in the range of 
interest) 

j<mg<fii<t: Aqcd < -, (12) 

where L is the spacial extent of the lattice and a is the lattice spacing. In this particular 
scenario, the scale Aqcd sets both the IR confinement scale that observables, such as nucle- 
ons and rho mass, are proportional to and the UV scale when the asymptotically free theory 
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becomes perturbative. As a result, these theories are thought to be phenomenologically 
similar to QCD. 

Another possible class of confining theories for a different set of parameter (often with a 
greater number of flavors), usually referred to as walking theories, follow a different hierarchy 

y < "^g < Ai/ < < Kuv < -• (13) 

Ij (X 

Most notably, the primary difference is the separation of the IR confinement scale and the 
UV perturbative scale. For a given renormalization scheme, the running of the coupling 
slows between the two scales, and this slow running is believed to give an enhancement to 
the chiral condensate proportional to {Kjjv / ^ir)" ^ where z/ is a power determined by the 
dynamics and renormalization scheme . For this reason, extracting the chiral condensate 
in the chiral limit (or, in terms of the xPT parameters, B/F) is the main priority for theories 
in the confining regime. 

For both classes of confining theories, xPT extrapolations are valid provided both and 
III are small enough. However, how small they have to be for a valid chiral extrapolation 
is expected to be significantly different between the two theories. There are two primary 
reasons for this. First, walking theories usually require a large number of flavors. Powers 
of the factor Nf appear in the coefficients of the higher order chiral logarithms and these 
coefficients can lead the higher order terms to be competitive with lower order terms if 
the pion mass is not small enough. This would result in a breakdown of the chiral EFT. 
Second, the pion mass itself is enhanced in these theories, so an even smaller quark mass is 
required to make a pion mass in the chiral regime. As a result, multiple lattice calculations 
of "confining" theories have found that chiral extrapolations do not hold for their current 
quark masses [60] (if running slow enough, such theories may resemble an IR conformal 



theory where anomalous dimension is fixed, 7(/i) ~ 7*, and the scaling analyses in Sec. IV 
may prove informative). It should be emphasized that a small isospin chemical potential 
(compared to the chiral breaking scale) is required for the following relations to hold. 

The study of the effects of isospin chemical potentials for QCD has been studied exten- 
sively for masses and condensates EH E2] and the Dirac spectrum This section 
will reflect a small subset of these effects as it pertains to extracting useful information for 
QCD-like theories. 



A. Pion mass 



The effect of the isospin chemical potential on the pion mass is separated in two regions, 
when the isospin chemcial potential is too small to form a pseudoscalar condensate {fi'j < 
2Bmq) and when it is large enough that a pseudoscalar condensate is formed (/Xj > 2Bmq). 
Both cases have been studied extensively EH E2] and the relevant tree-level results for 
extracting information about the condensate are presented here. Again, it should be noted 
that 2Bmq, fj,j <C A^ is required for the following analysis to be valid. 

In the flrst phase, fij < 2Bmq, no pseudoscalar condensate is formed and the chiral 
Lagrangian follows the form of Eq. ([2]). Here, much like zero density, the propagator matrix 
is has no cross terms (only terms proportional to ttq or 7r+7r_) and the particle masses are 
given by the zero values of its determinant 

detD-i = [p" -2Bmq][{j)- iiif -2Bmq][{j) + iiif -2Bmq] = Q, (14) 
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FIG. 2: Pion masses measured vs. (negative) isospin chemical potential. All axis are normalized to 
the LO vacuum pion mass value sJIBmq. The top curve (m^_), middle curve (mTrg), and bottom 
curve (?7i7r_^) are defined in Eq. (15) and Eq. (17). The curvature of after the pseudoscalar 
condensate forms (//| > 2Bmg) can give additional information useful for extracting the chiral 
condensate. 



and the resulting pion masses are (assuming the sign of fii is defined by fij = —\fii\) 

^2Bmq - , "^TTo = ^/'^Bruq , m^r. = ^/2Bmq + (15) 



In effect, at zero temperature, this phase only results in a shift in energy of the charged 
pions [HI]. While information about the condensate can be extracted from this phase, it is 
the same information that can be extracted from zero density simulations. 

In the second phase, fij > 2Bmq, a pseudoscalar condensate is formed and the chiral 
Lagrangian is now expanded about a new vacuum under the transformation S — )■ ^o^^o 
where So = ~ cos a + m ■ rsina. This leads to a propagator matrix that has no 
cross terms for the neutral particle and cross terms (proportional to 7r+7r+ or 7r_7r_) for the 
charged particles due to the flavor mixing of the pseudoscalar condensate. The zeros of the 
determinant (with three-momenta set to zero) are given by 

det = [pI - /i?] [pI - ^?(1 + 3(25m,/^2)2]p2 ^ (^g) 
and the resulting pion masses are given by 




, = , = l/i/U / 1 + 3 2~ • (!''') 



The first two masses do not contain any information about the chiral condensate, but the 
mass of the third pion is dynamically altered (as compared to the zero density case) as 
a function of the condensate and the isospin chemical potential. In Fig. [2| this can be 
interpreted as the curvature in when fij > 2Bmq. The usual Gell-Mann-Oakes-Renner 
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(GOR) relation {mlf^ = 2mq{%l)il))) no longer holds in the presence of the pion condensate, 
and the mass of m^r^ leads to a modified GOR-like relation given by 



f -y i2£2^(^"0' ^^^^ 

The primary advantage gained from this calculation (as pointed out in Ref.^Tj) is the 
addition of another "knob" besides niq to extract chiral behavior. Thus, for each within 
the chiral regime, multiple /i/ also within the chiral regime can be calculated (whose costs 
are determined by the symmetry breaking term e), and this information can lead to more 
control of chiral extrapolations. 



B. Pseudoscalar condensate 



In addition to extracting B/F from the pion mass, another approach would be to focus 
on measuring the pseudoscalar density directly. This has been performed stochastically 
in both quenched |13] and dynamical [Hl - liG] QCD at finite isospin density. In stochastic 
calculations of the chiral condensate, one nuisance is the presence of a large divergent piece 
on the order of the cut off scale, which masks the finite physics of interest. This can be seen 
in the naively discretized free propagator loop 



tr 



1 



_ I (27r)4 sin _ 



2m„ 



1 



2m„ 



_i (27r)4 sin^ap^ - 



(19) 



where the momentum dependence in the numerator is zero due to the trace over the 7-matrix. 
For the pseudoscalar condensate, a similar lattice artifact given by 



tr 



d^p 



(27r)^ ^7^ sin ap^ + niq + ie'y^r'^ 



2e 



d'^p 



1 



(27r)'^ ^ sin^ ap^ + + e 



2e 



(20) 



This divergent contribution will be smaller than its chiral condensate counterpart as long 
as e < niq. Neither divergence should depend on /i/ and consequently, taking differences at 
different /i/ may prove valuable. 



Following up Eq. (11), the pseudoscalar condensate at LO in the chiral expansion 
(m, /i/ ^ A^) is given by 



F^B sin a = F^B, 



( 2fim„ 



(21) 



The resulting pseudoscalar condensate as a function of isospin chemical potential for fixed 
niq and various values of B is shown in Fig. [3j Ultimately, we will want to take niq to zero 
(more specifically, take 2Bmq/ fij to zero). Upon doing this, the pseudoscalar condensate 
reduces to the quantities of interest, namely 



-r 



lim 75^^) 



F^B. 



(22) 
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FIG. 3: Pseudoscalar condensate (p.s.c) vs. isospin chemical potential (both normalized to a fixed 
small niq) for theories with B/rUq values of 1 (green), 3 (orange), 9 (magenta) and 27 (blue). As 
the chiral condensate increases, higher ^u/ values are required to form a pseudoscalar condensate. 

It should be noted that this particular result does rely on the fact that <^ A^, and 
corrections on the order of /xf are to be expected. However, as long as /x/ is well within the 
chiral regime, this particular method gives a clear method of extracting B from multiflavor 
lattice calculations with the large lattice artifact contributions potentially reduced. Also, 
getting at the chiral value for the pseudoscalar decay constant can be just as hard (if not 
harder) than extracting the chiral condensate. Thus, by combining this analysis with the 
pion mass analysis from the previous sections, B and F can be determined separately. 

IV. IR CONFORMAL THEORIES 

Upon reaching a certain critical value for the number of flavors in a theory with a given 
number of colors and fermions in a given representations, the theory is believed to be IR 
conformal due to the presence of an IR fixed point. While the theory still remains asymp- 
totically free for large energies, the fixed value for the coupling prevents any IR confinement 
scale from emerging. As a result, the theory does not have any spontaneous symmetry break- 
ing, and in the absence of any explicit massive scale, all observables have zero mass. For 
this reason, low-energy effective field theories that depend on spontaneous chiral symmetry 
breaking, such as xPT , are no longer valid. In the case of a mass deformed IR conformal 
theory, where a small quark mass is added (m^ ^ Ajyy), all hadronic masses follow the 



where 7* is the scheme-independent mass anomalous dimension. Many phenomenological 
arguments for understanding the quark flavor hierarchy from technicolor models that are 
"nearly" IR conformal (7(/x) ~ 7*) are based on the magnitude of 7* ^65| and for this reason, 
it is the primary quantity of interest for IR conformal theories. However, many intricacies 



relation [211 ESHS] 



Mo ~ m 1+^* , 



(23) 
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occur when extending this to the lattice. The first is that the relation only holds for small 
quark masses such that higher order terms are negligibly small. Another related concern is 
that in addition to the mass deformation, the inverse of the finite spacial extent, L~^, acts as 
an IR cutoff to the theory and the hierarchy L^^ <^ should be maintained to minimize 
the contributions from unphysical UV artifacts. 

Upon including the isospin chemical potential, the relevant hierarchy of scales in the 
conformal window are given by 

Y<mg<fii<^Auv< -■ (24) 

As will be shown in the subsequent sections, the implementation of the isospin chemical 
potential does not receive any renormalization due to the fact that it is coupled to a conserved 
current (the operator it is coupled to has engineering dimension as shown in Ref. [S3])- This 
feature allows for both varying the fermion mass near zero and allows for one to define finite 
density scaling relations akin to the finite size scaling. 



A. Universal scaling function 

The use of finite size scaling relations remains the best method for defining the universal 
scaling behavior and extracting the critical exponents for theories with an IR fixed point. 
This fact has not escaped the lattice community and the scaling behavior has been worked 
out for multiple classes of observables and operators [52] - [5l] . While finite size scaling has 
been explored on the lattice [23 EHl E31 ES], it ultimately proves to be quite demanding 
numerically since it requires multiple ensembles of different parameters while keeping 
small and <^ a~^. 

A similar set of universal scaling curves can be derived for the isospin chemical potential 
since the isospin chemical potential does not change under renormalization group (RG) 
transformations as a result of chemical potentials being coupled to conserved currents in the 
Lagrangian. The analysis for this "finite density" scaling (assuming the volume is large and 
fixed) is essentially the same as the analysis performed for finite size scaling in Ref. [53] , with 
the key difference being the spacial extent L is replaced with the isospin chemical potential 
fij^. Nevertheless, the key elements to this analysis will be repeated here. 

The correlation function calculated on the lattice is given by 

Coit;mg,fxj,L,ix) = (O^(t)O(O)) ^ Ae-^'^\ (25) 

where O is an artibrary interpolating operator and /i is the renormalization scale (not to be 
confused with the isospin chemical potential, ^i) and the r.h.s represents the long Euclidean 
time behavior of the correlation function. Assuming a fixed value for L with -C rriq, fij 
and defining rhq = nig/ /i, one arrives at the correlation function Co(t; rhq, fij, /i). Under the 
RG transformation fi ^ fi', where n = bn', the correlation function becomes 

Co{t;mg,fxi,fx) = r2^''Co(t; 6^+^*m„ /i,, /i')- (26) 

where 70 is the anomalous dimension of the field operator O. Performing a naive dimensional 
scaling by b to this result, 

Co{t;mg,fij,fi) = b-^^''^+'''^'>Co{b-H-b'+''*mg,bfij,fi). (27) 
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Now choosing b such that bj2j = Uq where Uq is a dimensionless constant. This leads to the 
correlation function 




Co{t■.m,,^^u^^)= \ — \ Co\ ^o, /ij , (28) 

where x = m^/ jj,]'^'^ and from the asymptotic behavior of the correlation function at long 
Euclidean time, the corresponding mass is given by a function of x, 

Mo ~ fiiF{x). (29) 

In the same way varying the box size is used in finite size scaling techniques, the relevant 
information of the mass anomalous dimension can be extracted by varying fij as the fermion 
mass approaches zero. In particular, when plotting Ma/fii vs. niq/ jj]^'^ at different /i/ 
values, the anomalous dimension can be determined by varying 7* and observing when all 
data points fall on the same curve. While this is akin to finite size scaling, the isospin 
chemical potential has more resolution at a comparable cost (can be adjusted ala quark 
mass without critical slow down) than the volume, which usually requires additional lattice 
sites (expensive, discrete) or tuning the lattice spacing (additional lattice artifacts). 



Using similar techniques, volume scaling can be incorporated into Eq. (29). Following 
the same procedure as above, any mass will now be a universal function of two variables, 
Mq ~ niF{x,y), where x = niq/ 11]^^ and y = hjL. 



B. Condensates 

As in confining scenarios, condensates can give important information regarding IR con- 
formal theories. In particular, the study of condensates in the presence of a chemical poten- 
tial in a given theory can provide another test as to whether or not a theory is conformal 
or confining. In an IR conformal theory with /x/ = and ruq = 0, all observables with 
mass dimension, including chiral and pseudoscalar condensates will be zero (technically, all 
observables should scale with since the finite spacial extent is a relevant operator in 
an IR conformal theory). The interesting question is the behavior of these observables at 
non-zero (but small) values of /x/ and rriq, and, in particular, how they compare to their 
confining counterparts. 

While spontaneous symmetry breaking does not occur in IR conformal theories, sponta- 
neous symmetry breaking also does not occur without a small source for confining theories 
at finite volume. Thus, as explained in Sec. [Tl| a small explicit symmetry breaking param- 
eter (which we called e) must be introduced to the system for the pseudoscalar condensate 
at fij > 2Bmq to be measured in both confining and IR conformal theories. However, the 
phenomena between these two theories should behave differently when the limit e — )■ is 



taken first. From Eq. (21), it is evident that as ruq — )■ 0, the pseudoscalar condensate is 
proportional to the confinement scale (in particular, the chiral condense). This is to be ex- 
pected, since spontaneous flavor symmetry breaking is expected to occur in infinite volume 
for theories with a large enough isospin chemical potential. However, the more interesting 
aspect stems from the fact that as rriq — 0, the pseudoscalar condensate has a diminishing 
dependence on /i/. 



12 



Within the conformal window, the pseudoscalar condensate takes a significantly differ- 
ent form when niq and /i/ are small, but non-zero, where m^^^^^ ■* plays the role of the 
confinement scale. Comparing the condensates in the chiral limit, 

hm .(V^^75^) ^ 1^'''; Confining Theory, _ ^^^^ 
mq^o 2 \Dij,j,IR Conformal Theory, 

where D is a dimensionless constant that could be zero. While this quantity is not an 
order parameter (unless D = 0), there is a difference in the behavior of the condensate as 
the isospin chemical potential is varied. By reducing the fermion mass to small values, the 
analysis of the pseudoscalar condensate is expected to display different behavior at several 
different (but small) isospin chemical potential values. Once again, it should be made clear 
how much this process depends on orders of limits. First, the explicit pseudoscalar breaking 
parameter e needs to be taken to zero, followed by niq, while keeping /i/ small enough such 
that these results will not be contaminated from higher order effects. 

V. ESTIMATES OF LATTICE ARTIFACTS 

In addition to the increased computational cost of inversions from reducing the quark 
mass, one always needs to be aware of the lattice artifacts that can arise. For example, 
simply decreasing the quark mass leads to rriT^L decreasing as well, which can introduce 
significant lattice artifacts, especially if m^L < 4. It should be noted that while this is 
a limitation for confining theories, large volume effects can act as a feature for finite-size 
scaling in IR conformal theories (more defined volume behavior). In addition to large IR 
effects, significant UV effects can become important, such as the emergence of non-physical 
phases such as the Aoki phase EB] • In this section we will address both of these examples 
when the quark mass is dropped in the presence of a small isospin chemical potential. The 
following two sections will address these lattice artifacts for the confining scenario only. 

A. Finite Volume 

The dominant volume effects in a lattice system with periodic boundary conditions are 
set by the mass of the lightest particle mass times the finite spatial extent of the system 
(for most theories, this is set by the product niT^L). In a confining QCD-like system, for 
small enough quark masses, xPT can be used to estimate the volume effects due to the finite 
spacial extent [69j. These effects emerge in the calculation of quantum corrections and loop 
diagrams, which the continuous momenta integrals are replaced by sums. In particular, for 
a periodic box, 

L 

The resulting finite volume effects for the pion propagator are often exponentially suppresses, 
as in the one loop correction term that appears in the calculation of the pion mass and chiral 
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where the asymptotic form of the Bessel function Ki{x) ~ e'^j^fx. As the number of 
flavors is increased, the coefficient of this exponential is also enhanced. In the presence of 
an isospin chemical potential, the volume effects are subtle. In particular, when \ii S> m^, 
two of the modes have ~ /i/, leading to heavier particles. However, the more noteworthy 
case is that in the condensed phase, one of the pseudoscalar particles becomes a massless 
Goldstone mode, whose correlation length is infinite and will always have volume effects. In 
practice, the mass of this Goldstone mode is cut off by the explicit flavor breaking term e 
that is deflned in Eq. (|9]) leading to ~ e. The corresponding volume effects depend upon 
the quantity eL, where there are two separate limits. The limit eL — )■ oo represents taking 
the box to inflnity before taking the explicit flavor breaking to zero. This limit is likely too 
computationally expensive to ensure. The second limit eL — )■ while all other mass scales 
times the box goes to inflnity will always be sensitive to the size of the system. However, 
if the box size is large enough, the system should be in the "epsilon regime" where the 
zero-mode volume effects are enhanced. Thus, epsilon regime chiral perturbation theory or 
random matrix theory can accurately account for these volume effects [TTHT^ . Additionally, 
mesonic two-point functions in the presence of an isospin chemical potential were previously 
calculated in the epsilon expansion [18] . 

B. Lattice spacing and Aoki regime 

The other well known lattice phenomena that occurs when reducing the mass while 
keeping all other parameters flxed is the emergence of the Aoki regime [671 EH] , where flavor 
and parity are spontaneously broken. This unphysical phase has been shown to depend on 
regions of the vriQ-g^ phase diagram and can be explained via a xPT argument based on 
estimates of the lattice spacing effects [HH]- This phase can play a signiflcant role when the 
LO xPT effect (given by the leading order pion mass) becomes on the same order as the NLO 
0{c?^ effect assuming the 0{a}^ effect comes in with the opposite sign. These competing 
effects, as worked out in Ref. [68], are most apparent when minimizing the potential of the 
xPT Lagrangian with the parameterization, E = A -|- iB ■ cr 



where C\ ~ and C2 ~ a^. In most lattice calculations with unphysical large pion masses 
and reasonably small lattice spacings, C\ 3> and these two effects never compete. However, 
as the pion mass is dropped while the lattice spacings are flxed, scenarios when c\ ~ C2 can 
occur and new unphysical minima can emerge in the potential. However, just because 
Ci ~ C2 does not mean that one is guaranteed to be in the Aoki regime when the actual 
lattice simulations are performed. 

When an isospin chemical potential is included, as ruq is dropped, rriT, ~ /i/ or ~ e. 
Before, we needed to ensure that our lattice spacing was chosen such that ^ o?hf^. 



V{A) = -ciA + C2A\ 



(33) 
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Now, /if ^ a^A^ and Ae ^ a^A^ is needed to ensure no Aoki phase emerges. It might 
appear that only fij ^ a^A^ is needed to have Ci 3> C2, and while this may end up being 
be the case, the determination of the sign and magnitude of Ci and C2 is a subtle non- 
perturbative question (i.e. large strange quark in lattice QCD does not necessarily prevent 
Aoki regime). It should also be noted that the Aoki regime and the condensed finite isospin 
density phase are quite similar in appearance, as they both produce a flavor and parity 
violating pseudoscalar condensate. Thus, lattice simulations in the Aoki regime may be 
incorrectly interpreted as an isospin condensed phase. 

Another aspect involving lattice spacing effects is the fact that irrelevant operators can 
lead to new mass-like terms in the Lagrangian proportional to aW where W ~ A^. Thus, 
since the LO combination in xPTis given by Brrig + aW, simply taking rriq to zero is not 
enough to remove the "lattice mass" . However, this lattice spacing effect should be reduced 
for improved actions, chiral lattice actions, or smaller lattice spacings. Nevertheless, negative 
values of niq may need to be used to acquire the "m^ — > 0" behavior discussed throughout 
this paper. 

VI. CONCLUSION 

In this work, we have presented a new scenario for extracting information from strongly 
coupled gauge theories by implementing an isospin chemical potential. Due to its positive 
definite determinant, this isospin chemical potential can be included in lattice caknilations 
and has multiple properties that can be used to help extract desired lattice quantities. First, 
when the isospin chemical potential is larger than its critical value, a pion condensate is 
formed. This pion condensate has multiple properties in the chiral regime that are beneficial 
for extracting the chiral condensate. Second, the isospin chemical potential couples to a 
conserved current and consequently, does not receive renormalization. To that end, a finite 
isospin density scaling analysis can be employed to extract the mass anomalous dimensions 
using similar techniques to finite size scaling. This scaling relation should hold for IR 
conformal and could also be applicable to slow running theories near conformality where the 
anomalous dimension is roughly constant, 7(/i) ~ 7*. Additionally, multiple ensembles with 
different chemical potentials can be simulated at comparable costs (in contrast to increasing 
volume), which gives improved resolution to both chiral or scaling analyses. These properties 
make a convincing case that this scenario should be pursued on the lattice. 

While the benefits of employing an isospin chemical potential arc plentiful, there are still 
multiple caveats that must be understood. The first, and most significant, caveat is the fact 
that much of the analysis in this work depends greatly on the fact that isospin chemical 
potential is in the chiral regime (/i/ <^ Am, Auv)- On the other end, for confining theories, 
volume effects are now governed by the combination /iiL and eL, and analytical techniques 
in the eL — )■ limit might be needed. Also, contributions from Aoki phases are also possible 
and should be taken into account. 

One could make the observation that the inclusion of an isospin chemical potential while 
taking rUg and e to zero is simply moving the previous issues with the quark mass to the 
chemical potential. In many ways, this astute point is indeed correct and for that reason, 
this method should be viewed as another tool for probing lattice phenomena. However, the 
two other key features of the isospin chemical potential, namely the formation of a pion 
condensate and coupling to a conserved current, make this system inherently different from 
the usual lattice system with just a quark mass. The majority of the interesting physics for 
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these systems not only occur in the chiral hmit, but greatly depend on how the theory varies 
near the chiral limit. With the presence of the isospin chemical potential, these systems can 
be explored in exactly this way. 
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